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Abstract 



p. 

f-n | Experiment shows that the reverse length-lexicographical word ordering 

consistently yields far smaller Grobner bases for modular p-group algebras 
than the length-lexicographical ordering. For the so-called Jennings word 
ordering, based on a special power-conjugate group presentation, the as- 
sociated monomial algebra is a group invariant. The package Present finds 
■ Grobner bases for these three orderings. 
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These notes document some experiments I did in the late 1990s to compare the 
usefulness of local and well orderings for noncommutative Grobner basis calcula- 
tions in modular group algebras of p-groups. They are not (or rather no longer) 
intended for publication. They are provided in the hope that the tables in sec- 
tion [5] may be of interest. 

Part of the material here formed the basis for Chapter 1 of my book0. The 
system Present is no longer being maintained, but parts of it live on in the p-group 
cohomology packaged for Sage. 
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1 Introduction 



Let G be a finite p-group. One way to do computational homological algebra 
over the modular group algebra ¥ P G is to realise it as a finitely presented asso- 
ciative algebra, together with a Grobner basis for the relations ideal. There 
are then several methods available for computing minimal projective resolu- 
tions jFeustel et al. 1993j . |E. Green et al. 1998] . |D. Green 1999] . 

The algebra generators are required to lie in the Jacobson radical. Usually 
one chooses generators Qi for the group, and then takes the a, = gi — 1 as gener- 
ators for F p G. See |Farkas et al. 19 93J or [D. Green 1997] for an introduction to 
noncommutative Grobner bases. 

The choice of word ordering is of prime importance in the design of efficient 
Grobner basis methods. This paper compares the usefulness in our context of 
three word orderings, and announces a package Present for computing Grobner 
bases with respect to these orderings. They are: 

• The length-lexicographical ordering (LL): the standard ordering for ho- 
mological algebra over associative algebras. 

• The reverse length-lexicographical ordering (RLL) allows one to read 
off the radical layers. Empirically, it always yields a small Grobner basis. 

• The Jennings ordering (see Definition 13. 5j) is defined on special presen- 
tations, which are constructed using the Jennings series and involve some 
redundant generators. Broadly speaking, it is for p-group algebras what 
power-conjugate presentations are for p-groups. All Grobner bases have 
the same size, and in fact the associated monomial algebra is an invariant 
of the group's order. 

The criteria for comparing orderings are: 

• Smallest size of Grobner basis, taken over all choices of generators. 

• How easy is it to realise a small Grobner basis? That is, how sensitive is 
size of Grobner basis to choice of generators? 

• How many properties of the Grobner basis (e.g., maximum length of a 
reduced word) do not depend on choice of generators? 

In Sections [2] and [3] we look at the orderings one by one. For each ordering, 
a method is presented for computing the Grobner basis for the relations ideal. 
These methods were implemented in the package Present. Data structures and 
other aspects of the implementation are discussed in Section HI Results obtained 
using Present are presented in Section [5], where we draw conclusions about the 
comparative usefulness of the orderings. 

The package Present is written partly in Magma code and partly in C. See 
the end of the paper for details of how to obtain it. 
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2 Grobner bases 

First we recall the basics about Grobner bases in free associative algebras. Let M 
be the free associative monoid with 1 on a finite set A. That is, the elements 
of M are words in the alphabet A. Assume we have an ordering on M which is 
admissible in the following sense (far weaker than that of [Feustel et al. 1993] ): 

Definition 2.1 Let M be a free monoid. An ordering on M is called admissible 
if uviiu < uv 2 w whenever u, Vi,v 2 , w are elements of M with V\ < v 2 . 

Now let A be a set of algebra generators for ¥ P G, all lying in the Jacobson 
radical J(W P G). Then M is the set of monomials in the free associative F p - 
algebra on A. Write / for the relations ideal: the kernel of projection from this 
free algebra to ¥ p G. A word w G M is called a tip if it is tha largest word in the 
support of some element of /; if not, it is called a nontip. That is, the nontips are 
the words whose images in ¥ P G are linearly independent of the images of their 
predecessors. It follows immediately that the images of the nontips are linearly 
independent in ¥ P G. 

Lemma 2.2 With the above assumptions, the nontips form a basis for ¥ P G. 

Proof. By Nakayama's Lemma, J N (¥ p G) is zero for some N. All algebra 
generators lie in J, and so only finitely many words are nonzero in ¥ p G. These 
span ¥ p G, and contain the nontips as a spanning subset. ■ 

By admissibility, the tips form an ideal in the free monoid M. This ideal has a 
unique smallest generating set. It consists of the minimal tips, a tip being called 
minimal if and only if all proper subwords are nontips. Define 

Q := {w — p(w) | w a minimal tip} , 

where u(w) is the linear combination of nontips which in ¥ p G equals w. Then 
Q and the length iV words together generate I. So, in the presence of the con- 
straint A N C J, the set Q is a Grobner basis for /: the unique completely reduced 
Grobner basis for this ordering. 

Remark 2.3 For the length-lexicographical and Jennings orderings, the con- 
straint A N C / is not necessary. 
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2.1 The length-lexicographical ordering 



Denote by £(w) the length of a word w G M. Putting an ordering on the set A 
of algebra generators induces a lexicographical ordering <i cx on M. The length- 
lexicographical ordering < L l on M is then defined as follows: 

W\ <ll w 2 if i{wi) < £{w 2 ), 

or same length, and W\ <i cx w 2 . 

This ordering is admissible. Moreover, each word w G M has only finitely many 
predecessors under < LL . We do not use the Buchberger algorithm to compute 
the Grobner basis, because for p-group algebras it takes an unreasonable amount 
of time to stop. Rather, we exploit the fact that ¥ P G is finite-dimensional, which 
means that it is feasible to list all the nontips. 

Computing a Grobner basis Construct G in Magma using a faithful per- 
mutation representation. Choose minimal generators g 1 ,...,g r for G, and set 
= 9i — 1- Take the at as the set A of algebra generators. Record the matrix 
for the right multiplication action of each Oj, with basis the group elements. 

Using these matrices, each word is constructed as a linear combination of 
the group elements, starting with the smallest word and proceeding in LL-order. 
Gaussian elimination allows us to distinguish the tips from the nontips. We stop 
when we have found all the nontips: their number is known in advance. The list 
of minimal tips is deduced from the list of nontips by word manipulation. Change 
of basis gives us the matrix for the multiplication action of each generator with 
respect to the basis of nontips, and we can read off v{w) for each minimal tip w 
from these matrices. 

2.2 The reverse length-lexicographical ordering 

Again defined on the free monoid M on an ordered set A, this admissible ordering 
is the opposite of LL. Namely, 

Wi <rll w 2 if and only if w 1 > LL w 2 . 

Note that RLL is not a well-ordering. 

Lemma 2.4 The KLL-nontips of length at least r constitute a basis for J r (¥ p G). 
Hence the length r RLL-nontips are a basis for a complement of J r+1 in J r . 

Proof. The words of length at least r span J r . If a linear combination of 
length r words lies in J r+1 , then the largest length r word involved is a tip. ■ 

Magma can compute the smallest iV such that J N is zero (using the Jennings 
series, the dimensions of all radical layers can be computed). So we could obtain 
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the nontips by running through the words of length at most N in RLL order. But 
this could take forever. For example, if G is a Sylow 2-subgroup of the sporadic 
finite simple group C03, then N = 23 and A has size 4. Hence the number of 
words of length N is 2 46 . 

Computing a Grobner basis A length r word is a tip if and only if it lies 
in the space spanned by its length r RLL-predecessors and by J r+1 (W P G). Using 
Gaussian elimination and the matrices for the algebra generators, we can calculate 
a basis for J r+1 from one for J r . So we start with r = and work through to 
r = N. All length r nontips are words of the form w.a with a G A and w a length 
r — 1 nontip. We work through these words in RLL-order. Once we have the 
nontips, we proceed as for LL. 

3 The Jennings ordering 

For r > 1, define F r (G) to be the rth dimension subgroup of the finite p-group G. 
That is, 

F r (G) = {geG\g-leJ r (¥ p G)}. 

Clearly G = Fi(G) > F 2 {G) > • ■ ■ . The F r do eventually reach the trivial group, 
and they form a central series, the Jennings series for G. This series need not be 
strictly decreasing, however. These and more facts about the Jennings series are 
demonstrated in Section 3.14 of [Benson 1991] 

Definition 3.1 Let G be a p-group of order p n . Elements g±, . . . ,g n of G are 
Jennings pc-generators for G if the first mi elements are minimal generators 
for Fi(G)/F 2 (G), the next m 2 are minimal generators for F 2 /F 3 , and so on. 

Jennings pc-generators do yield a polycyclic presentation for G, as we now see. 

Lemma 3.2 Let gx, . . . , g n be Jennings pc-generators for G. For each 1 < i < n, 
we have gf G (g i+1 , . . . , g n ); and [g u gj] G (g j+1 . . . , g n ) for 1 < i < j < n. 

Proof. It is known that [F r , F s ] < F r+S and that pth powers of elements of F r 
lie in F pr . m 

Example 3.3 Let G = (a, b, c, <fi) be the following semidirect product group of 
order 32: the subgroup (a, b, c) is normal, and elementary abelian of order eight. 
The order of is four, and the effect of conjugation on the left by <fi is 

a 1 — > b 1 — > c 1 — > abc . 

Then F 2 = (ab, ac, 4> 2 ) and F 3 = (ac) are elementary abelian, and so a, 4>, ab, 4> 2 , ac 
are Jennings pc-generators for G. 
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The last three generators are Jennings pc-generators for F 2 . Another family 
of Jennings pc-generators for F 2 is ac,ab,<p 2 . This family cannot be extended 
to Jennings pc-generators for G, since [<ft, ab] = ac. Hence to find Jennings pc- 
generators for G, it is not enough to take minimal generators for G and add on 
Jennings pc-generators for F 2 . 

Definition 3.4 Let gi, . . . ,g n be Jennings pc-generators for G. Then ai, . . . , a n 
are Jennings generators for the group algebra ¥ P G, where — gi — 1. We say 
has dimension r if aj G J r (¥ p G) \ J r+1 (¥ p G), or equivalently gi G F r \ F r+ \. 

Order the Jennings generators as follows: a\ < a 2 < ■ ■ ■ < a n . For a word w in 
these generators, define its dimension dim(u>) in the obvious additive way: 



Definition 3.5 The Jennings ordering is the following ordering on the monoid 
of words in the Jennings generators: 



Observe that the Jennings ordering is admissible; and that dim(if) > £(w) for 
each word w. The next result tells us that it is straightforward to calculate the 
Grobner basis for the Jennings ordering. 

Proposition 3.6 Let CLi , . . . , CL n be Jennings generators for ¥ P G. The mini- 
mal tips are the words a? { and with j < k. The nontips are the words 
affa^Si ■ ■ ■ a!] 1 with < ej < p — 1. A nontip w lies in J r \ J r+1 if and only 
dim(w) = r. 

Proof. We just have to show that a\ and a^a^ (j < k) are tips. For then all 
nontips are of the form claimed, and the number of nontips equals the number 
of claimed nontips. So the nontips and hence the minimal tips are as claimed. 
The distribution in radical layers then follows from Jennings' theorem (Theorem 
3.14.6 in [Benson 1991J V 

If dim(aj) = r then a\ lies in ¥ p F pr , which is generated by the aj with 
dim(aj) > pr. Hence a\ is a linear combination of words in these aj. Each 
such word either has dimension greater than pr or has length one. So af is a tip, 
since it has length p. 

Similarly, if dim(aj) = s and dim(afc) = t, then c = [gj,gk] lies in F s+t . Set 
7 = c — 1, so that 7 G ¥ p F s+t . Arguing as above, 7 is a linear combination of 
words smaller than ajat or a^Oy Since g^g^ = cg^gj, we have 



r 




W\ <j w 2 if dim(wx) > dim(w 2 ), 

or dimensions equal and £(wi) < £(w 2 ), 

or same dimension and length, and w\ >i ex w 2 . 



aja^ = a^Oj + 7 + 70^ + 70, + ^a^j 
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and the right hand side is a linear combination of words smaller than a^. ■ 

Example 3.7 Let G be the cyclic group of order four. Then Jennings pc- 
generators for G are gi,g 2 with g\ = g 2 , g\ = 1 and [gi,g 2 ] = 1- The cor- 
responding Jennings generators for ¥ p G are a± = g± — 1 and a 2 = gi — 1. The 
minimal tips are af , aia 2 and a|; the nontips are a 2 ai, a 2 , a 1; 1 in ascending order; 
and the Grobner basis is 

a\ + a 2 , a x a 2 + a 2 a x , a\ . 

Remark 3.8 Jennings pc-generators are usually not minimal group generators. 
This means that there are relations involving length one words. However, the 
Jennings ordering does guarantee that all tips have length at least two. 

Computing a Grobner basis Get Magma to compute the Jennings series. 
Use this to pick Jennings pc-generators for G. The nontips and the minimal tips 
are known by Proposition 13.61 Proceed as for LL. 

4 The implementation 

In this section we provide an overview of the package Present, and describe the 
data structures used. Some components of Present are written in Magma code, 
others are written in C. The C components use M. Ringe's C MeatAxe to handle 
vectors and matrices over finite fields. 

Groups must be constructed in Magma as permutation groups. The function 
regularPermutationAction is provided to convert pc-groups (and hence matrix 
groups too) into permutation groups using the regular permutation action. 

4.1 Selecting minimal generators 

To choose minimal generators for a p-group G in pc-presentation, Present first 
asks Magma for the Frattini subgroup $(GQ, and sets H = Then it selects 

an element g G G — H, adds this to the list of generators, and replaces H by 
(g,H). This is repeated until H is equal to G. At this point the elements g 
constitute a minimal generating set. 

There remains the question of which element of G — H to select at each stage. 
Two generator selection methods were investigated: 

• The most obvious method: pick g G G — H completely at random. Minimal 
generators constructed in this way are called arbitrary minimal generators. 

• Calculate the exponent of each element of G — H, and pick g at random 
from amongst those with the smallest exponent. Such minimal generators 
are called smallest exponent minimal generators. 
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The merits of these two methods are compared in Section [5J By default, Present 
uses the smallest exponent method for LL, and the arbitrary method for RLL. 

4.2 The nontips tree 

The group generators are passed to the C programs as a list of permutations. For 
the Jennings ordering, generator dimensions are passed too. The Grobner basis 
is now determined by the method for the chosen ordering given above. 

The nontips are words, and all proper subwords of a nontip are again nontips. 
Hence the most natural way to store the nontips is as a tree, in which the children 
of a word w are the nontips of the form w.a with a e A, the set of algebra 
generators. On the other hand, the nontips constitute an ordered basis, and so 
we also want to be able to address them as an array. 

In Present, the nontips are stored as a length \G\ array of pathnodes. The 
class pathnode contains 

• the word w itself, its array index and its length 

• an array of pointers to its children, indexed by the elements of A. (If w.a 
is a tip, then the corresponding pointer is null.) 

• a pointer to the parent word. Also, which child of the parent it is. 

The nontips are stored in the array in ascending order for LL, and in descending 
order for RLL and the Jennings ordering. This means that the root of the tree 
is always located in the first entry of the array; and that the nontips tree can be 
built while the nontips are being determined. 

For the second point, observe that the LL nontips are determined in ascending 
order. Also, the number of length r RLL nontips is known in advance, and they 
are determined in ascending order after all nontips of smaller length have been 
determined. Hence for LL and RLL, the index of each nontip is known the 
moment it is identified as a nontip. 

Strictly, the Jennings nontips are an exception here. They are written down 
in lexicographical order, and then sorted into Jennings order. 

4.3 The components of Present 

The package Present consists of the package MakeBasis of Magma functions, and 
four C programs. 

The Magma function makeBasis constructs nontips and action matrices for 
the LL and RLL orderings. A variant has a user-defined number of attempts at 
finding the smallest Grobner basis. The function makeJenningsBasis constructs 
nontips and action matrices for the Jennings ordering. There is no need to have 
repeated attempts. 
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For each ordering, a flag allows the user to insist that the defining group 
generators are used. This eliminates the random component. 

The C programs makeNontips and makeActionMatrices are invoked by the 
Magma functions to determine the nontips and the action matrices respectively. 
The program writeGroebnerBasis can then be used to write out the Grobner 
basis and the minimal tips. The utility program grouplnf o prints out relevant 
statistics on the groups by decoding the nontips file header. 

5 Results and conclusions 

The package Present was used to compare the orderings LL and RLL. The groups 
used in the comparison were: all 51 groups of order 32; all 267 groups of order 64; 
and eight Sylow p-subgroups of sporadic finite simple groups. 

Define LL(G) to be the smallest size of a Grobner basis for the relations ideal 
with respect to the LL ordering. Define RLL(G) similarly for RLL. Using Present 
we can obtain empirical approximations eLL(G) and eRLL(G) to these numbers. 

By Propostion 13.61 all Grobner bases for the Jennings ordering have the same 
size Je(Cr). This is \n{n + 1) for a group of order p n . 

5.1 Groups of order 32 

There are 51 groups of order 32. The smallest Grobner basis found in twenty 
attempts was recorded for each combination of: group, ordering (LL or RLL) 
and generator selection method (arbitrary or smallest). 

For the RLL ordering, the generator selection method made no difference. 
For LL, smallest minimal generators yielded a smaller Grobner basis in three 
cases: one less element for two groups, and three less for one group. 

The comparative performance of the two orderings is shown below. The most 
extreme difference was for the group with Hall-Senior number 43, where eLL = 31 
and eRLL = 10. 



d := eLL — eRLL 


d < 


d = 


1< d < 3 


4 < d < 6 


7 < d < 9 


d > 10 


No. of groups 





17 


16 


8 


7 


3 



5.2 Groups of order 64 

There are 267 groups of order 64. The smallest Grobner basis found in twenty 
attempts was recorded for each combination of group, ordering and generator 
selection method. 

For the RLL ordering, smallest minimal generators yielded a larger Grob- 
ner basis in 63 cases, and a smaller Grobner basis in one case. Each time, the 
difference was only one element. For LL, smallest minimal generators yielded a 
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Arbitrary 


Smallest exponent 


Group 


Order 


Min 


Max 


n 


cr 


Min 


Max 


n 


a 


Sy\ 2 (HS) 


2 9 


128 


338 


261 


50 


104 


127 


114 


5 


Syl 2 (M 24 ) 


2 io 


433 


906 


706 


101 


183 


543 


373 


81 


Syl 2 (Co 3 ) 


2 io 


502 


871 


690 


92 


212 


574 


405 


85 


Syl 3 (McL) 


3 6 


126 


333 


258 


50 


133 


245 


197 


43 



Table 1: Generator selection methods compared for LL. Each sample size: 30. 





Arbitrary 


Smallest exponent 


Group 


Order 


Min 


Max 


li 


a 


Min 


Max 


li 


a 


Syl 2 (HS) 


2 9 


10 


13 


11.0 


1.1 


11 


13 


11.9 


0.9 


Syl 2 (M 24 ) 


2 io 


15 


22 


17.3 


1.5 


15 


20 


17.1 


1.0 


Syl 2 (Co 3 ) 


2 io 


13 


21 


15.5 


1.7 


13 


21 


16.2 


2.0 


Syl 3 (McL) 


3 6 


11 


14 


11.8 


1.0 


12 


14 


12.6 


0.9 



Table 2: Generator selection methods compared for RLL. Each sample size: 30. 

smaller Grobner basis in 97 cases, with mean difference 3.3; and a larger Grobner 
basis in 69 cases, with mean difference 4.6. 

The comparative performance of the two orderings is shown below. The most 
extreme difference was for the group with Hall-Senior number 187, where eLL = 
57 and eRLL = 10. 



d := eLL — eRLL 


d < 


d= 


1 < d < 8 


9 < d < 16 


17 < d < 32 


d > 33 


No. of groups 





29 


112 


73 


50 


3 



5.3 Sylow ^-subgroups of sporadic finite simple groups 

One of the main applications driving the development of group cohomology soft- 
ware is calculating the cohomology rings of sporadic finite simple groups. So 
Sylow p-subgroups of these groups are important examples for Present. More- 
over, they are a good source of larger p-groups. 

The sensitivity of Grobner basis size to generator choice was investigated 
by looking at four such Sylow p-subgroups. The program was run thirty times 
for each combination of: group, ordering (LL or RLL) and generator selection 
method. The results are shown in Tables [1] and [2j We give the smallest and 
largest sizes of Grobner basis found, together with mean and standard deviation. 

In Table El we compare smallest Grobner basis sizes for all three orderings 
by looking at eight Sylow p-subgroups of sporadic finite simple groups. Each 
empirical value is based on at least twenty calculations. 
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Sylow 2-subgroup of 




Syl 


3 Of 




M 22 


HS 


M 24 


Co 3 


Suz 




McL 


Suz 


Order 


2 7 


2 9 


2 io 


2 io 


2 13 


2 14 


3 6 


3 7 


eLL 


34 


104 


150 


236 


2669 


3111 


126 


417 


eRLL 


8 


9 


15 


13 


21 


18 


11 


13 


Je 


28 


45 


55 


55 


91 


105 


21 


28 



Table 3: Smallest known Grobner bases. Each sample size at least 20. 
5.4 Conclusions 

Grobner bases for the length-lexicographical ordering LL are consistently much 
larger than for the reverse length- lexicographical ordering RLL. In Table EJ the 
size of the smallest LL Grobner basis lies between 14% and 33% of the group 
order. This means that the LL ordering is unsuitable for p-group algebras. 

By contrast, RLL consistently yields very small Grobner bases: the size be- 
having very roughly as the logarithm of the group order. Moreover, RLL allows 
one to read off the radical layers in the group algebra. Hence the RLL ordering 
is well-suited for computing with p-group algebras. 

The Jennings ordering also seems to be suitable for computing with p-group 
algebras. Finding Jennings generators is easy, and it would appear that all Jen- 
nings generating sets are equally useful. Many properties of a Grobner basis are 
invariants of the group's order, including the size. Grobner bases are not as small 
as for RLL, but this is compensated for by the fact that fewer multiplication 
operations are necessary to obtain the action of the average minimal tip. 

The smallest exponent generator selection method seems to deliver a small LL 
Grobner basis more often than the arbitrary method does: see Table [U However, 
there are numerous exceptions amongst the groups of order 64. For the RLL 
ordering, all evidence suggests that the smallest exponent method should be 
avoided: it is no more reliable at yielding a small Grobner basis than the arbitrary 
method, and sometimes misses the smallest Grobner bases. 
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